NASA-TM 

74091 Turbulent Bounda/y Layers Over 
NonsTationary P,llane Boundaries 

y 

NASA Technical Memorandum 74091 


Turbulent Boundary Layers Over 
Nonstationary Plane Boundaries 


Alan T. Roper and Garl L. Gentry, Jr. 
FEBRUARY 1978 




NASA Technical Memorandum 14091 


Turbulent Boundary Layers Over 
Nonstationary Plane Boundaries 


Alan T. Roper 

Rose-Hulman Institute of Technology 
Terre Haute, Indiana 

and 

Garl L. Gentry, Jr. 

Langley Research Center 
Hampton, Virginia 


IVIASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
information Office 


1978 



SUMMARY 


Methods of predicting integral parameters and skin-friction coefficients 
of turbulent boundary layers developing over moving ground planes are evalu- 
ated using test information from the Langley V/STOL tunnel. These data were 
obtained at Reynolds numbers higher than those previously available for evalu- 
ating these predictive techniques. The three methods evaluated are (1) rela- 
tive integral parameter method; (2) relative power law method; and (3) modified 
law of the wall method. 

Methods (1) and (2) can be used to predict moving-ground-plane shape fac- 
tors with an expected accuracy of ±10 percent. They may also be used to predict 
moving-ground-plane displacement and momentum thicknesses with lower expected 
accuracy. This decrease in accuracy can be traced to the failure of approxi- 
mations upon which these methods are based to prove universal when compared with 
V/STOL tunnel test results. 

Although no attempt is made to establish the accuracy of the local skin- 
friction coefficients predicted by the various methods, those produced by 
methods (1) and (2) show a high degree of agreement. The modified law of the 
wall method predicts a more rapid decrease in skin friction with increasing 
ratio of ground-plane velocity to free-stream velocity than do methods (1) 
and ( 2 ) . 


INTRODUCTION 

Turbulent boundary-layer flows developing over nonstationary boundaries 
are of interest in studies related to a variety of flows including the Ludwieg 
tube, moving shock waves near solid boundaries, and the near field drag of tube 
vehicles. It is thus desirable to accurately predict the boundary-layer char- 
acteristics and the skin-friction coefficient for such flows. 

Since all present methods of turbulent-boundary-layer calculation rely upon 
empirical information, and since no significant body of such information exists 
for the moving-ground-plane case (except for refs. 1 to 3), any useful method 
of prediction must originate from existing stationary-ground-plane methods 
(ref. 4). Three such methods, (1) relative integral parameter method, (2) rela- 
tive power law method, and (3) modified law of the wall method, have been pro- 
posed in references 1 to 3 and partially verified by tests performed in two 
Langley Research Center wind tunnels equipped with moving ground belts. The 
bulk of this test information was obtained at Reynolds numbers less than 2 x 10^. 
To extend the Reynolds number range and complete the verification of one or more 
of the proposed methods, the present tests were conducted in the Langley V/STOL 
tunnel at Reynolds numbers up to 9.2 x 10^. In this report results of the 
latest series of tests are compared with those predicted by means of the three 
proposed methods. Where appropriate, results of the lower Reynolds number tests 
(refs. 1 to 3) are also included in the comparison. 



SYMBOLS 


The units used for physical quantities in this report are given in both 
the International System of Units (SI) and the U.S. Customary Units except that 
stations on figures are given in Units only. Measurements and calculations 
were made in the U.S. Customary Units. 

A,B,A',B' empirically determined constants (eqs. (15) and (18)) 

Cf total skin-friction coefficient 

Cf(0) stationary-ground-plane total skin-friction coefficient 

Cp local skin-friction coefficient 

Cf<(0) stationary-ground-plane local skin- friction coefficient 

g(R) empirically determined polynomial (eq. (6)) 

H shape factor, <S*/9 

H(0) stationary-ground-plane shape factor 

N reciprocal of power in relative power law formulation (eq. (8)) 

N(0) reciprocal of power in stationary-ground-plane power law formulation 

n reciprocal of power in power law formulation (eq. (7)) 

R velocity ratio (nominal value), Vg/U 

Rjj Reynolds number based on x and U 

Ry Reynolds number based on y and U 

U x-component of mean velocity at edge of boundary layer, 

m/sec (ft/sec) 

u local x-component of mean velocity, m/sec (ft/sec) 

u^ friction velocity, ^T^/p, m/sec (ft/sec) 

Vg velocity of ground belt, m/sec (ft/sec) 

X longitudinal distance from leading edge of ground belt, m (ft) 

y vertical distance from gound belt, cm (in.) 

6 boundary-layer thickness (defined as y at which u/U = 0 . 99 ), 

cm (in.) 
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DISCUSSION 
Experimental Data 

The data presented in this report were obtained in the Langley V/STOL 
tunnel. The tunnel was fitted with a moving smooth-surface ground belt (rep- 
resenting a moving ground plane) and operated with the tunnel side walls and 
ceiling in place. Slots fitted in the side walls to allow variable porosity 
were, however, partially opened. The general arrangement of the test facility 
is shown in figure 1(a). Velocity profiles were determined from total-pressure 
measurements obtained with a total-pressure rake composed of 46 separate tubes 
of 0.1 27- cm (0.05 -in.) diameter. Simultaneous static-pressure measurements 
were made by means of a separate probe located to one side of the rake and at 
a nominal height of 4.2 cm (1.65 in.) above the ground belt. The general 
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arrangement of the total pressure rake is shown in figures 1(b) and 1(c). All 
pressure data were obtained by means of a rotary valve and pressure transducer 
arrangement . 

The V/STOL tunnel is equipped with a porous suction plate to remove the 
normal tunnel floor boundary layer ahead of the belt; therefore, boundary layers 
presented in this report are assumed to originate at the leading edge of the 
belt. In practice, however, it is not possible to remove the entire tunnel 
boundary layer. The effects of this incomplete boundary-layer removal are shown 
for the R = 1.0 (nominal) profiles taken in the V/STOL facility in figure 2. 
Additional disturbances to the flow were caused at the leading edge by the small 
slot (between the suction plate and the belt) and the natural entrainment of 
air from this slot by the belt motion. Inaccuracies were encountered for read- 
ings taken nearest the 'belt because of the tendency of the belt to lift slightly 
when in motion. None of the data are corrected for these effects. 

Total-pressure measurements were made at a nominal dynamic pressure of 
526.7 Pa (11 Ibf/ft^), velocity of 29.3 m/sec (96.00 ft/sec), and Reynolds num- 
ber per unit length of 1.8? x 10^ per m (0.57 x 10^ per ft) at stations 1.37 m 
(4.50 ft), 2.62 m ( 8.58 ft), 3.35 m (11.00 ft), 4.27 m (14.00 ft), and 5.18 m 
( 17.00 ft) behind the belt leading edge. Tunnel static pressure during testing 
was nominally atmospheric. For these conditions, the range of test Reynolds 
numbers was 2.62 x 10^ ^ R^ = 9.21 x 10^. Data were recorded at values of 
Vg/U of approximately 0, 0.24, 0.48, 0.74, and 1.00. 

Throughout the test a periodic unsteadiness in the boundary layer of as 
much as ±0.8 m/sec (±2.5 ft/sec) was noted in the velocity measurements. No 
such fluctuations were apparent in the high-velocity central core of the 
test section where tunnel-center-line velocity varied less than ±0.2 m/sec 
(± 0.50 ft/sec) during any given test. Constant monitoring of ground belt 
speed failed to reveal fluctuations of a frequency sufficient to explain the 
periodic unsteadiness observed in the flow. All efforts to find and eliminate 
the source of the disturbance were unsuccessful. To eliminate errors induced 
by this flow condition, multiple readings were taken of all pressure measure- 
ments as summarized in table I. Velocities quoted in this report are the arith- 
metic averages of the velocities determined for each of the multiple pressure 
measurements. This scheme was not totally satisfactory and some of the scatter 
apparent in the measured velocity profiles is undoubtedly attributable to this 
source. 

Velocity profiles measured during the V/STOL tunnel tests are presented 
in figure 3; the velocity data are also given in table II. Summaries of the 
boundary-layer thickness 6 and the integral parameters 6*, 9, and H 
determined from the measured data are presented in table III. The variation 
of the boundary-layer thickness with Vg/U is presented in figure 4. Also 
displayed for comparison in that figure are results from the lower Reynolds 
number tests reported in references 1 and 2. 
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Theoretical Methods 


Brief descriptions are presented of the three methods suggested (refs. 1 
and 2) for predicting the development of the turbulent boundary layer over a 
moving ground plane. Evaluation of each method is based on a comparison with 
data measured during the V/STOL tunnel tests and, where appropriate, with data 
from references 1 and 2. 

Relative integral parameter method .- This method (see refs. 1 and 2) is 
based^upon the development of expressions for the relative integral parameters, 
6*, 9, and H, obtained by replacing the variable u/U in the usual parameter 

definitions (ref. 4) by 


u u - Vb (u/U) - R 

0 ■ U - Vb " 1 - R 


( 1 ) 


Thus, for example 



If the following assumptions (1) H is nearly independent of R and 
(2) 0/0(0) is some function independent of x, say g(R) are adopted, 

expressions relating to stationary and moving-ground-plane integral param- 
eters can be generated 


0 = g(R) 9(0) 


(2) 


H 


H(0) 

1 + R[H(0) - l] 


(3) 


H(0) 

6* = g(R) 0(0) (4) 

1 + R[H(0) - lj 

where g(R) is obtained by fitting the measured data. 

The accuracy of this method rests ultimately upon the validity of assump- 
tions (1) and (2). Of these two assumptions, the first is the least critical. 
Errors in predicted values of H/H(0) resulting from deviations in the assumed 

relation H/H(0) = 1 can be expressed as 
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e(H/H(0)) = ± 


(5) 


1 - R 

[l + R(H - 1)]2 


e(H/H(0)) 


The variation of H/H(0) as a function of R is presented in figure 5 for the 
three sets of test data . It will be noted that assumption ( 1 ) is rather seri- 
ously violated (e(H/H(0)) maximum of the order of ±0.25). The most serious 
variations correspond to large values of R. However, as the first term of 
equation (5) is always less than one, it serves to attenuate errors resulting 
from assumption (1). Further, the magnitude of that attenuation grows with 
increasing values of R. For data in this report, maximum errors in H/H(0) 
resulting from the assumption H = H(0) are no more than ±10 percent.- These 
errors correspond to midrange values of R. Errors in the determination of 
H/H(0) also enter into the determination of 6*/6*(0). (See eq. (4).) 


The accuracy of assumption (2) is more critical as it enters directly 
(unattenuated) into both the determination of 0/0(0) and 6*/6*(0). (See 
eqs. (2) and (4).) Reference 2 suggested that g(R) could be represented 
adequately by the expression 


g(R) = (1 - R)[l - 1.1756R + 0.7863R(1 + R)] (6) 

The variation predicted by equation (6) is displayed in figure -6 together 
with data from the V/STOL tunnel tests and references 1 and 2. The expres- 
sion appears to be a reasonable representation of the data, although a greater 
x-dependence is apparent in the current data than in that of either of the 
previous tests. Scatter about the predicted variation indicates a spread of 
approximately ±0.10 which is nearly uniform across the R range. Some of this 
scatter is undoubtedly due to the unsteadiness noted in the boundary-layer 
velocities of the V/STOL tests. 

Measured values of H/H(0), 0/0(0), and 6*/6*(0) are compared with 

those predicted from equations (2), (3)> and (4) in figure 7. Comparisons 
between measured and predicted values of H/H(0) are good, the maximum error 
being less than ±10 percent. Comparisons for 6*/6*(0) and 6/9(0) are poor 
and refle’ct the scatter apparent in g(R) (fig. 6). 

The additional assumption that Cf/Cf(0) (like Cf/Cf(0)) is indepen- 
dent of X allows the local skin-friction coefficients for the stationary and 
moving-ground-plane flows to be related. 

Relative power law method .- The power law method (see ref. 1) depicts the 
velocity within the boundary layer as 



where n is a parameter somewhat dependent upon Reynolds number. Examination 
of existing moving-ground-plane boundary-layer measurements (refs. 1 and 2) 
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led Roper and Gentry (refs. 1 and 3) to conclude that velocity distributions 
measured relative to a moving boundary admit to a similar formulation, 



where N is nearly independent of Vg/U. Although data from the current series 
of tests display more scatter than those of earlier tests (ref. 1), they conform 
sufficiently to the format suggested by equation (8) to corroborate its validity 
at least far downstream of the start of the belt (fig. 8). 

Slopes of all measured profiles (V/STOL tunnel tests, refs. 1 and 2) are 
summarized in figure 9 in the forms suggested by both equations (7) and (8). 

It is important to note that although n is strongly dependent upon Vg/U, 

N is more nearly constant. The primary motivation for the power law formula- 
tion is its simplicity. To maintain that simplicity, it appears desirable to 
set N = N(0) and to accept the resulting inaccuracy implied by the Vg/U 
dependence apparent in figure 9. 

The usual integral parameters can be developed by using their respective 
definitions and equation (8) 


6 » = 


(1 



(9) 


6 


(1 



e(o) 


( 10 ) 


H 


N + 2R) 


H(0) 


( 11 ) 


Once stationary-ground-plane integral parameters have been computed, equa- 
tions (9), (10), and (11) can be used to predict the corresponding moving- 
ground-plane values provided a suitable expression for 6/6(0) can be 
developed. Measured values for 6/6(0) are presented in figure 10. Much 
scatter is apparent, particularly at large values of Vg/U. The scatter is 
undoubtedly due partially to the inherent difficulty in determining 6 and 
6 accurately, and, for the V/STOL series of tests, is aggravated by the 
unsteadiness encountered in the boundary layer. 

Measured values of 6/6(0) and the values of N(0) obtained from least- 
squares fits of the stationary-ground-plane logarithmic profiles were used with 
equations (9), (10), and (11) to compute 6*/6*(0), 9/0(0), and H/H(0). 

These predicted values are compared with the experimentally determined ratios 
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in figure 1 1 . Comparisons did not appear to be impressive but were acceptable 
and display increasing accuracy of prediction as the longitudinal distance from 
the belt leading edge increased. 

For incompressible turbulent flow past a semi-infinite flat plate, Blasius 
(ref. 4) gives the variation of Tq with 6 as 


To 

= 0.0225 


pu2 



( 12 ) 


Extending equation (12) to include relative velocity (ref. 5), replacing 1/4 
by 1/n, and rearranging yield the expression for local skin-friction coef- 
ficient (ref. 1) 


2n-1 

Cf Tq / y \1/ri 

— = = 0.0225(1 - R) n / -—1 (13) 

2 pu2 ypU6j 

where RSI. Equation (13) can be used to predict the local skin-friction 
coefficient for moving-ground-plane flows once the parameter q has been 
determined. Substitution of equation (13) into the integral form of the momen- 
tum equation for a flat-plate zero-pressure-gradient boundary layer, integra- 
tion, and rearrangement of the results produce the expression 


n 

f 0.0225 n + i\n+i 

6 = 1 (1 

\ 0/6 n / 


1 

2n-i — 



(14) 


By using measured values of 6 and 9/6 and appropriate values of the various 
test conditions, equation (14) can be solved numerically for the values of q. 
These results are presented in figure 12 for both the present V/STOL tunnel 
tests and those of reference 1. The variation of q/q(0) from both series of 
tests is consistent. Equation (14) was solved for- Cf by using computed values 
of q. The resulting local skin-friction coefficients are presented in fig- 
ure 13 together with results from reference 1. With the exception of a single 
data point from the latest tests, scatter in the predicted value of cp as a 
function of Vg/U is less than ±13 percent. 


Modified law of the wall method .- The final approach suggested in refer- 
ence 1 is based upon the hypothesis that the relative velocity profiles obey 
the law of the wall in the modified form 


u 



+ B 


(15) 
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where A = 5.75, B = 5.20, and 


ut 




Cf 

T 


( 16 ) 


The local skin-friction coefficient based upon absolute velocity is related 
to that based upon relative velocity by the equation 


Cf Gf 

_ = (1 - r)2 _ (17) 

2 2 


Adopting the formulation suggested by Clauser (ref. 6), equation (15) can be 
rewritten as 

; = A* log (R„) + B' (18) 

U 

where A’ and B' are functions of A, B, and Cf. Charts constructed for 
a range of Cf values using equation (17) and measured data from the V/STOL 
tunnel tests are presented in figure 14. Comparison of measured data with the 
computed Cf grid is used to determine the value of relative skin-friction 
coefficient. 

Although some of the profiles lack a clearly defined linear portion, it 
was nonetheless possible to determine approximate values of cp for all sta- 
tions except X = 2.62 m (8.58 ft). Values of Cf determined in the manner 
just outlined are presented in figure 15. Also displayed in that figure are 
the results reported in reference 1 . Agreement between results for the two sets 
of measured data is good, maximum scatter being less than ±6 percent. 


, Comparison of Local Skin-Friction Coefficients 

The relative integral parameter method, relative power law, and modified 
law of the wall methods all yield means for predicting local skin-friction 
coefficient. The local skin- friction coefficient is an extremely difficult 
quantity to measure directly over a moving ground plane and no attempt has been 
made to establish the accuracy of the various methods. Results from both the 
V/STOL tunnel tests and those presented in reference 1 are displayed in fig- 
ure 16. Note that predictions of the relative integral parameter and power 
law methods agree within less than ±6 percent. Values of cp predicted by 
the modified law of the wall method are generally lower than those of the 
other methods. 
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Comparison of Predicted and Measured Integral Parameters 


In this section comparisons are made of measured values of 6*/6»(0), ' 
0/6(0), and H/H(0) with those predicted by the relative integral parameter 
and power law methods. Where available, data from the lower Reynolds number 
tests as well as from the V/STOL tunnel tests are displayed. 

Shape factor parameter (H/H(0)) .- Expressions for H/H(0) were developed 
in the relative integral parameter (eq. (3,)) and relative power law methods 
(eq. (11)). The first expression contains the approximations that H = H(0), 
whereas the second approximates the actual value of N by N(0). 

The accuracy of equations (3) and (11) is displayed in figure 17 by plot- 
ting predicted against measured values of H/H(0). Lines indicating ±10 percent 
error are included in figure 17- The relative integral parameter method is the 
somewhat more accurate of the two methods yielding errors from -9 percent to 
+2 percent. This method generally somewhat underestimates the actual value of 
the ratio H/H(0). With the exception of four measurements from the current 
set of tests, the relative power law yields errors from -3 percent to +13 per- 
cent. The power law method generally overestimates values of H/H(0). 

Momentum thickness parameter (6/6(0)) .- Expressions for 6/0(0) were 
developed in the relative integral parameter (eqs. (2) and (6)) and relative 
power law methods (eq. (10)). Measured values of 6/6(0) were taken from the 
curve fitted to the measured data in figure 10. 

Once again, the accuracy of equations (2) and (10) is depicted by plotting 
measured against predicted values of the ratio 6/6(0) (fig. 18). Lines 
representing ±10 percent and ±20 percent errors are included in that figure. 

Of the two methods, the relative integral parameter method appears most accu- 
rate, 80 percent of the predicted points involving errors of ±20 percent or 
less and nearly 50 percent of the predictions falling within ±10 percent of the 
measured value. For the relative power law method, 56 percent of the predicted 
values are within ±20 percent and 38 percent within ±10 percent of the measured 
values . 

The accuracy of the proposed methods in the determination of 6/0(0) dete- 
riorates markedly from that displayed in the prediction of H/H(0) (fig. 17). 
This deterioration is due in large part to the failure of the empirical expres- 
sions for g(R) and 6/6(0) to accurately represent the experimental data. For 
the relative power law predictions, the situation is further aggravated by the 
approximation that N = N(0). 

Displacement thickness parameter (6*/6*(0)) .- Expressions for 6»/6*(0) 
were developed in the relative integral parameter (eq. (4)) and relative power 
law methods (eq. (9)) where, once again, values of 5/6(0) are taken from a 
curve fitted to the measured values of that' parameter (fig. 10). 

The accuracy of equations (4) and (9) is depicted by plotting measured 
against predicted values of the ratio 6*/6*(0) (fig. 19). Lines representing 
±10 percent and ±20 percent errors are included in that figure. Each of the 
two equations predict 60 percent of the measured data points within errors of 
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less than or equal to ±20 percent. The relative power law method predicts 34 
percent of the data points within errors of less than or equal to ±10 percent 
whereas the integral parameter method predicts only 26 percent within this 
error band. Overall, there is less dispersion of the error chart data than that 
displayed in the 0/0(0) chart (fig. 18). Once again, the lack of an accurate 
match between measured data and the empirical formulations for g(R) and 6/6(0) 
seems to be the major source of error. 


CONCLUSIONS 

The validity of the proposed methods is somewhat clouded by the scatter 
apparent primarily in the Langley V/STOL tunnel test results presented in this 
report. This scatter is believed to result from the unsteadiness noted in the 
boundary-layer velocity measurements for which no explanation is advanced. 
Nevertheless, the data support the conclusions of NASA TN D-6788, at least far 
from the leading edge of the moving ground belt. These conclusions are; 

1 . The velocity profiles measured near a moving ground plane can be 
described by a relative power law formulation of the form 



where u/U is the relative velocity ratio, 6 is relative boundary -layer 
thickness, y is the vertical distance from the ground belt, and N (recip- 
rocal of power) is invariant with ground-plane speed. 

2. The integral parameters (displacement thickness 6*, momentum thickness 
0, total skin-friction coefficient Cf, and shape factor H) can be described 
with reasonable accuracy by either the relative integral parameter or the 
relative power law methods. 

3. The modified law of the wall method predicts a more rapid decrease in 
local skin-friction coefficient with increase in belt velocity to free-stream 
velocity ratio than do either the relative integral parameters or relative 
power law methods. 

As the result of the Langley V/STOL tunnel tests an additional conclusion 
can be drawn. 

4. Although no statement can be made concerning their accuracy, there 
is a high degree of agreement between the variations of skin-friction ratio 
(Cf/Cf(0)) with momentum thickness ratio (6/9(0)) predicted by the relative 
integral and relative power law methods. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 10, 1977 
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TABLE I.- SUMMARY OF MULTIPLE DATA POINT MEASUREMENTS 


X 

R = Vg/U for - 

ro 

ft 

10 data points 

30 data points 

1.37 

2.62 

3.35 

4.27 

.18 
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8.58 
1 1 .00 

14.00 

17.00 

None 

0, 0.243, 0.487 
0, 0.242, - 
0, 0.242, 0.486 
-, 0.242, 0.484 

0, 0.245, 0.491, 0.739 

-, -, -, 0.737 

-, -, 0.492, 0.740 
None 

0, -, 0.685 




TABLE II.- DATA FROM BOUNDARY-LAYER SURVEYS 
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TABLE II.- Continued 
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(c) General arrangement of total pressure rake. AH dimensions 

in centimeters. 




Figure 2.- Velocity profiles displaying effect of incomplete removal of normal 
wind-tunnel floor boundary layer ahead of moving ground belt. 
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Figure 3-- Measured velocity profiles. 
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Figure 3.- Continued. 
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Figure 3.- Continued. 
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Figure 3-- Concluded. 



Figure 4.- Variation of boundary-layer thickness with ground-plane speed. 






Figure 6.- Variation of g(R) = 6/0(0) with ground-plane speed. 
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Figure 7.- Continued. 
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Figure 7.- Concluded. 
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Figure 8.- Measured relative velocity profiles. 








.242 



42 


Figure 8.- Concluded. 
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(a)- Station x = 1.37 m (4.50 ft). 

Figure 11.- Comparison of measured and predicted values of integral 

relative power law method. 
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Variation of reciprocal of exponent in Blasius shear stress expression as 
a function of ground-plane speed. 
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Figure 14.- Skin- friction coefficients determined by modified law of wall method. 
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Figure 14.- Continued. 
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Figure 1^.- Continued. 
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Figure 14.- Continued. 
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Figure Concluded. 
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Figure 16.- Comparison of local skin-friction coefficients determined by- 
relative integral parameter method, relative power law method, and 
modified law of wall method. 
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H(R)/H(0), experimental 






6*(R)/6*(0), experimental 





